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ABSTRACT 

We study the properties of satellites in the environment of massive star-forming galaxies at z ~ 1.8 in the COSMOS field, using a 
sample of 215 galaxies on the main sequence of star formation with an average mass of ~ 10** Mq. At z > 1.5, these galaxies typically 
trace halos of mass > lO*** Mq. We use optical-near-infrared photometry to estimate stellar masses and star formation rates (SFR) of 
centrals and satellites down to ~ 6 x 10* Mq. We stack data around 215 central galaxies to statistically detect their satellite halos, 
finding an average of 3 galaxies in excess of the background density. We fit the radial profiles of satellites with simple /i-models, 
and compare their integrated properties to model predictions. We find that the total stellar mass of satellites amounts to ~ 68% of the 
central galaxy, while SED modeling and far-infrared photometry consistently show their total SFR to be 25 - 35% of the central’s rate. 
We also see significant variation in the specific SFR of satellites within the halo with, in particular, a sharp decrease at < 100 kpc. 
After considering different potential explanations, we conclude that this is likely an environmental signature of the hot inner halo. 
This effect can be explained in the first order by a simple free-fall scenario, suggesting that these low-mass environments can shut 
down star formation in satellites on relatively short timescales of ~ 0.3 Gyr. 

Key words. Galaxies:halos - Galaxies:evolution - Galaxies:high-redshift - Galaxies:star formation 


1. Introduction 


Although the gradual infall of small dark matter halos onto 
larger ones has become a relatively straightforward aspect of 
the standard hierarchical formation paradigm, what happens to 
the baryons they contain is less well understood. In particular, 
the mechanisms that drive the evolution of their constituent 
galaxies become more complex as they are accreted by larger 
structures. Of special relevance are the processes that regulate 
and ultimately suppress star formation in galaxies in the early 
Universe. Their relationship to, and influence on, the galaxies’ 
immediate environment is not known with certainty; also 
debated is the relative importance of internal mechanisms versus 
externally driven ones (although the former are expected to 
be dominant in massive systems and the lat ter to act preferen¬ 
tially on lower-mass satellite galax ies; e.g.. iBaldrv et alJbOOfil : 
iPeng et aTll2010l: iGabor et^uOlIh . The z = 1.5 - 2.5 epoch 
is particularly interesting as a transition period when global 
star formation in the universe peaks, but also where the first 
ostensibly collapsed and virialized galaxy structures appear, 
in which a spatial segregation of different galaxy types (e.g., 
passive and active) is observed. In particular, the cores of mas¬ 
sive clusters appear to bec ome dominat e d by quiescent galaxies 
around this time (e.g., ISnitler et al.l 120121: IStrazzullo et al.l 


1201 31 : lOobat et al.ll2013h . From a theoretical point of view, the 
increasing temperature of the gaseous medium in group- and 
cluster-scale halo s starts to efficiently prev e nt accretion aroun d 
this epoch (e.g., iDekel & Birnboirnl 120061 : iDekel et ^ l2009h . 
thus affecting the build-up and evolution of the galaxies they 
host. One can therefore expect the processes regulating mass 
accretion onto galaxies, crucial to our understanding of galaxy 
build-up, to be relatively accessible to observation at this epoch. 


For practical and historical reasons, the mass regime most 
often explored at high redshift has been that of large galaxy 
clusters, which are the richest and most readily selectable halos. 
They are also the most biased regions in which to study environ¬ 
mental effects on galaxy evolution. However, at high redshift the 
advantages offered to galaxy evolution studies by their galaxy 
density and mass contrast are somewhat counterbalanced by 
their relative rarity. A lot of effort has thus been devoted to the 
search for high-redshift structures, through a variety of methods. 
Although remarkable progress has been made recently, only 


a handful of z > 
terized so far (e.e 

1.5 structures have been accurately charac- 

.. Andreon et a 

1. 20091: iPanovichet all 120Id: 

iGobat et al .11201 it 

Stanford et al.l 

20121). as this endeavor is still 


fundamentally hampered by the relatively limited area for which 
deep datasets are available (although this may change thanks to 
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ongoing and future deep wide-field infrared surveys). On the 
other hand, the lower “group” mass range, at this redshift that of 
the progenitors of z = 0 clusters, has been less systematically 
explo r ed; it requires either much deeper data (e.g.. lErfaniar et alJ 
120131; iTanaka et al.l 12013h or tracers, such as specific galaxy 
types (e.g., quasars or giant radio-galaxies), that correlate with 
structures but are not directly proportional to mass (unlike, e.g., 
galaxy density or diffuse X-ray emission). This latter type of 
tracer has been historically used for higher r edshift systems, 
such as proto-clusters (e.g.. iHatch et al.l[201 it IWvlezalek et al.l 
l2013h . At lower halo masses, even isolated massive galaxies are 
expected to be at the center of galaxy assemblages, as a simple 
consequence of the hierarchical nature of matter distribution in 
the universe. 


In iBethermin et al.l (l2014li . we indeed found that massive 
(~ 10" Mq), star-for ming galaxies on the main sequence 

of st a r formation (e.g ., Brinchmann et al.l l2004t iDaddi el'iD 
l2007t iRodighiero et aijl2010l) in the range z - 1.5 - 2.5 had 
clustering properties consistent with halos of mass (> 10'^ Mq). 
Subsequent stacking of deep X-ray datasets available in COS¬ 
MOS yielded constraints on the total mass consistent with the 
clustering analysis (~ 2 x 10'^ Mq), thus providing independent 
confirmation. This suggests that massive main-sequence galax¬ 
ies constitute a conspicuous tracer of group-scale environments 
at z ~ 2 (more so than, e.g., isolated quiescent galaxies), thus 
allowing for the easy study of these systems. 


Small halos comprising a central galaxy and its satellite 
system are particularly useful for probing the environmental 
dependency of galaxy properties over large mass and redshift 
ranges. In particular, they provide a powerful tool to constrain 
quenching mechanisms and timescales, thro ugh their tell-tale 
signa t ure on mass profiles a nd functions (e.g..lWang et al.ll2010[ 
12014 iPhillips et al.l I20l4 [Hartley et al.l 1201^ or simple de- 
rived quantities suc h as the fraction of quiescent satellites (e.g., 
iGeorge et al.l201 iT) . Being vastly more abundant and structurally 
simpler than massive galaxy clusters, these systems allow for 
a straight forward test for galaxy assembly and evolution mod¬ 
els (e.g., IGuo et af] l2Ql ih . without requiring deep knowledge 
of their components (age, position in phase-space, etc.). Here 
we have taken an intermediate approach, focusing on the prop¬ 
erties of star-forming satellites, and in particular on the varia¬ 
tion of their star formation rates (SFR). We have only consid¬ 
ered systems with massive star-forming centr als; while quies¬ 
cent g alaxy pairs also trace similar sized halos dBethermin et al.l 
|2014|) , their center of mass is less clear, which would make an 
investigation of the radial dependency of satellite properties less 
straightforward. This paper is organized as follows; in Section 
|2] we describe the dataset and our sample selection. In Sec¬ 
tion |3l we present the integrated properties of satellites and the 
method used to derive them. In Section 0] we discuss their vari¬ 
ation with environment and present our conclusions in Section 
|5] Throughout this paper we assume a ACDM cosmology with 
Ho - 70 km_s“* Mpc“\ Q,„ = 0.27, and A = 0.73, and a 
IChabrieil (l2003h initial mass function (IMF; relations used here 
that assume a different IMF have been converted to this one). 
Magnitudes are given in the AB photometric system throughout. 


2. Data and sample selection 

In this work, we have used aperture-corrected photometry from 
the X,-selected catalog of the COSMOS/U ltraVISTA survey 
dMcCracken et al.l l20l2i from iMuzzin et al.l (l2013h . We have 



M. (Mo) 

Fig. 1. SFR of centrals as a function of their stellar mass, for SFR es¬ 
timates from, respectively, HerschelfPACS fluxes (red circles) and SED 
modeling (black squares). Our adopted parametrization of the main se¬ 
quence at z = 1.8 is shown as a gray line, with the 0.5 dex limit shown 
by dashed lines. The SFRs shown here have been corrected according 
to the difference in normalization between the main sequence at z = 1.8 
and the redshift of the objects. 


adopted photometric redshifts (Zphot) from lllbert et ^ (120131 
and references ther ein) rather than the Zphot estimates from 
iMuzzin et al.l (1201 3h . as the former had access to a larger and 
deeper training se t of spectroscopic redshifts (Zspec) (especially 
zCOSMOS DeeD; lFillv et'^l2()07l) . ensuring greater reliability 
of Zphot at z > 1. For objects for which th ese were not available , 
we have used the Zphot estimates from IMuzzin et al.l (l2013l) . 
Where possible, we have a l so us ed Zspec from the zCOSMOS 
Bright sample (iFillv et al.l I 2 OOI . Stellar masses, SFR, and 
rest-frame colors were then recomputed based on this merged 
catalog, as described in Section 12.11 Finally, we have also 
used mid- and far-infrared (FIR) maps from Spitzer/MlPS 


(iLe Floc’h et al. 

2009h. HerschellFACS (from the PEP survev: 

iFutz et al.l 2011 

), and Herschel/SPIRE (from the HerMFS 

survev; Oliver et nTl2()T^ in the analvsis. although onlv the 

PACS data were 
below). 

used for the construction of the sample (see 


For consistency with IBethermin et al.l (I20l4) . we have built 
a sample of mas sive, star-forming galaxies by considering 
all BzK-selected (iDaddi et al.l 12004 . i/eric/zel/PACS-detected 
objects with SF Rs within 0.5 dex of the main sequence (as 
parametrized in IBethermin et al.l 1201 2h . using the recomputed 
stellar mass derived from the UV-near-infrared (NIR) photom- 
etry and SFR derive d from fits to the 100 and 160jU fluxes with 
iMagdis et alJ (1201 2h templates. Selecting only PACS-detected 
sources mostly yields massive galaxies, while the 0.5 dex 
removes lower and upper outliers (e.g., quenching galaxies and 
starbursts, respectively). The position of galaxies in our sample, 
relative to the main sequence at z = 1.8, is shown in Fig. [T] 
In addition, we have also rejected galaxies that are less than 
1' (~ 500 kpc) away from known overden sities with redshifts 
consistent within the 6 8% confidence level dChiang et al.ll20l4 
IStrazzullo et alJ 120151) . or probable companions (also with 
redshifts consistent at 68% confidence) of mass mfrac ^ 2 
times higher within r,,,> = 35" (~ 300 kpc). This constraint 
corresponds roughly to the typical expected virial radius 
(r,,,>) of a ~ 2 X 10'^ Mq halo in the sample’s redshift range 
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(1-5 < z < 2.5), while the constraint on the mass ratio between 
the central and companion galaxy reflects the typical stellar 
mass uncertainty w hen all variables and degeneracie s are taken 
into account (e.g., iBerta et ^ 120041 : IConrovI l2013h . The first 
criterion is conservative and meant to minimize the risk of 
including halos significantly more massive than ~ 2 x 10'^ M©. 
Such systems would also be richer and might somewhat bias 
the results of our analysis. These values, although physically 
motivated, were also chosen as a compromise to obtain a 
relatively clean but still statistically significant sample. We note 
that altering them slightly (e.g., nifrac = 1 and r = 5" so as to 
reject probable interacting pairs) does not change the sample 
much, nor the outcome of the analysis presented below. This 
criterion yields 215 galaxies, which we henceforth assume are 
central to their halo (hereafter, centrals), with a mean mass of 
(M*) = 1.3 X 10" M© and a mean redshift of <z) = 1.8. The 
scatter of the centrals’ stellar mass distribution is ~ 0.25 dex, 
consistent wi th that expected from a single halo population at 
this redshift dBehroozi et alJl2()T^ . The distribution of masses 
and redshifts in the sample is shown in Fig.|2 


2.1. Stellar population modeling 


We have estimated stellar masses, SFR, and dust extinction for 
both centrals and satellites (see SectionO by fitting the available 
UV-N IR SEDs with two different types of iBruzual & Chari^ 
(l2003li models: stellar masses were computed assuming a 
generalistic delayed exponential star formation history (SFH) of 
the form SFR(t) = SFRq x (f/r^) x exp(-f/T) which includes 
both rising and declining cases. The age t and timescale t were 
allowed to vary between 0.1 Gyr (f) or 1 Myr (r) and the age 
of the Universe at the redshift of the galaxy. Star formation 
and extinction values, on the other hand, were estimated by 
fitting the rest-frame UV photometry assuming a constant SFR 
with an age limit of f > 100 Myr ( the timescale ex plicitly 
assumed for UV-derived SFRs; e.g., iKennicutl Il998h . This 
was done because, while the rest-frame optical-NIR part of the 
SED reflects the entire SEH of the galaxy, the rest-frame UV 
part is sensitive to “instant” SER. This also makes comparison 
with the literatur e easier, since direct UV-to-SER conversions 
(lKennicutl]ll998t) generally assume continuous star formation 
over ~ 100 Myr (see also Section]?]). We have included extinc- 
ti on by dust, consider ing values of E{B-V) < 2 and assuming 
a iCalzetti e t al.l (l20 00li functional form wi th an additional UV 
bump ( Noll et alT 20091: iBuat et all l2012h . This limit should 
be safely above the normal v alues for even the most massive 
centrals in our sample (e .g., iGarn & Best! l2010t IZahid et ^ 
l2014t IPannella et al.l]20f4]) . Solar metallicity was assumed for 
all models: due to the age-metallicity degeneracy, this parameter 
would only matter in the case of very old populations, which 
are not liable to be relevant for the galaxies we consider in 
our analysis. We have ignored the far- and near-UV GALEX 
bands, as the resolution of GALEX is significantly poorer 
than that of the other instruments contributing to the catalog 
(> 4" compared to sub-arcsecond seeing), which precludes 
efficient deblending on scales typical of the size of halo cores 
(~ 5", and thus relevant to the analysis in Section]?]). Eor the 
same reason, we did not include the 5.8 and 8 pm SpitzerflRAC 
bands in the modeling and treated the EIR data separately, as 
detailed in Section [TTI The SED modeling was thus performed 
on a maximum of 25 photometric bands. 


3. Stacking analysis and integrated properties 


Eor each central galaxy, we have constructed a sample of 
candidate members of that galaxy’s halo (hereafter “satellites”) 
by selecting all uncontaminated, non-stellar sources in the 
m erged catalog with X , < 23.4, the completeness limit cited 
by iMuzzin et al.] ( |2013b , and with zl 6S < Zcen < Zh6S- Here Zcen 
is the photometric redshift of the central and zun (zh6&) the 
lower (respectively upper) 68% confidence limit to the Zphot 
of the putative satellites. The uncertainty on the association 
between galaxies is likely to be dominated by the uncertainty 
on the photometric redshifts of the fainter satellites, rather than 
of the bright centrals. Accordingly, we assume the redshift of 
the centrals to be fixed at the best-fit value. In addition, we have 
only considered sources within a radius of 80", or ~ 700 kpc, 
about twice the typical r^ir of the host halos at this redshift. Since 
the average minimum separation between centrals in our sample 
is ~ 2', and the redshift range in which satellites are selected 
can be large, this has allowed us to estimate the contribution 
of background and foreground interlopers while minimizing 
cross-contamination between satellite systems. These satellite 
subsamples were then decomposed into passive and star-forming 
galaxies based on their rest-frame U-V and V-J colors (e.g., 
IWuvts et al.l l2007l hereafter UVJ). The mass distribution of 
satellites selected as star-forming and passive, respectively, is 
shown in Eig. ]2] The distribution of passive-selected satellites 
is shown for completeness only, as we focus on star-forming 
systems from Section]3T]on. 


Satellites were grouped in concentric annuli of width 
2" (~ 17 kpc) centered on each halo galaxy, and their rele¬ 
vant properties (photometric fluxes, stellar masses and SER) 
averaged in each radial bin. These values were then divided 
by the total area of the annulus. In each case, the contribution 
of background and foreground interlopers (hereafter, “back¬ 
ground”, for convenience) was estimated from the values in 
annuli between 50" and 80". This region was chosen to be 
comfortably distant from any significant galaxy excess due to 
the host halo (as seen in Eig. ]2) and cover a large area. We note 
that this method yields values consistent with background levels 
determined through random apertures. However, the use of an 
annular region that is still relatively close to the central should 
better account for local background variations due to interloper 
clustering (e.g., galaxy filaments). We use the same background 
radii for all subsamples and do not attempt to, e.g., adapt the 
size of the bins to the mass of the central, since r,„> is not very 
sensitive to halo mass variations and the stellar mass scatter of 
the centrals is consistent with that of a single halo population. 
Uncertainties in stacked properties were estimated, in each bin, 
using 1000 bootstrap resamplings of the data, with sizes of half 
the initial sample. However, bootstrap-derived uncertainties 
tend to be underestimated as they do not account for systematic 
uncertainties. This is typically more noticeable in the case of 
large samples. In an attempt to compensate for it, we rescaled 
the bootstrap estimates so that, when fitting the background as 
a constant term, the reduced chi-square dt Xq - 1 if initially 
larger. This is an ad hoc correction meant to produce more 
conservative error estimates (however, systematic errors are not 
necessarily Gaussian and the use of a estimate might not be 
formally justified). As an example, Eig. ]?] shows the averaged 
stellar mass, SER, B-, and Xj-band flux density profiles of 
satellites. 
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Fig. 2. Distribution of stellar masses (left) and redshifts (middle) of centrals in our sample, and of stellar masses of satellites (right) color-selected 
as star-forming (blue histogram) and quiescent (red histogram). 



Fig. 3. Left: Positions of all satellite candidates (open circles) compared to their respective centrals (filled circle), with the overdensity most visible 
at < 10". The size of each symbol varies as a function of its stellar mass. The red dashed contours show scales of 10", 20", and 30", respectively. 
Right: fraction of expected “real associations” among satellite candidates, as a function of radius (where tiscair) = n{r) — iii,, n(r) and n* 

being, respectively, the number density of satellite candidates at radius r and at > 50"). 


The average number of satellites per central, and the 
contribution of satellites to the total stellar mass and SFR 
of the halos, were estimated by fitting a parametric function 
to t he profiles. We have consider ed both the NFW pro- 
file (iNavarro. Frenk & White Il995h and projected j8-model 
dCavaliere & Fusco-Femian^ 19781) . commonly used for, re¬ 
spectively, dark matter halos and galaxy clusters. We find that 
the latter provides significantly better agreement (AAIC~ 45) 
with the data, especially at small (r < 5") radii where the NFW 
profile is too shallow. This might be the result of mild mass 
segregation, as shown in Section IS (see also, e.g ., Watson et al.l 
l2()12t Ivan der Burg et al.l 120141 : ll^cionere et al.ll2015 ). On the 


other hand, the profiles are well fit with /? ~ 0.9. Although a 
precise constraint of the true shape of the satellite profile is 
beyond the scope of this paper, this is co nsistent with values 
derived at both low an d high redshift (e.g.. lPoDesso et aDl2004l: 
IStrazzullo et al.lU013l) . As a check, we have also performed fits 
letting the background vary freely, which yielded background 
values consistent with those estimated from the outer annuli. 


Integrating the number density profile, we find the average 
excess of Ks < 23.4 satellites t o be 3.3 + 0.2. This is some¬ 
what above the value reported by iHartlev et al.l (1201 5l) for high- 
redshift centrals and could reflect the different nature of our sam¬ 
ple as well as the higher average mass of its centrals. Similarly, 
we find the integrated stellar mass and SFR of satellites to be 
M* = (3.6 + 1) X 10'° Mo and SFR = 28 + 6 Mq yr^', or re¬ 
spectively 28% and 15% of the average central mass and SFR 


(1.3 X 10" Mq and 192 Mq yr“*, respectively). However, these 
values are underestimates for two reasons: first, when selecting 
satellite candidates, we have considered objects with photomet¬ 
ric redshifts consistent with the central’s at only the 68% con¬ 
fidence level. This was done to minimize background contami¬ 
nation when estimating radial profiles and trends (see above and 
Section|4|i. However, assuming to the first order that the width of 
the intrinsic velocity distribution of satellites is negligible com¬ 
pared to photometric errors, and that the background redshift dis¬ 
tribution is mostly flat in the redshift range of the selection, we 
can expect to lose 32% of satellites to redshift uncertainties. The 
actual integrated mass and SFR should then be higher by a factor 
of 1.47. Secondly, we have only considered satellites within the 
completeness limit of the catalog, Ks - 23.4, and thus do not in¬ 
clude fainter satellites. At z ~ 2, the Ks band does not trace stel¬ 
lar mass equally for all galaxies, due to the flux contribution from 
young stars becoming non-negligible. Through comparison with 
stellar population models, we have estimated the corresponding 
mass limits for quiescent and star-forming galaxies to be, respec- 
tiv ely, logMvf=10.3 and 9.8. The models used here were based 
on iBruzual & Chariot (l2003l) templates assuming, respectively, 
a single burst of maximal age and a main-sequence SFH of the 
form 

( IVl \~0.2 

(l+z(0)' Meyr-' (1) 

where M(t) and z(t) are, respectively, the stellar mass and 
redshift at time t after the onset of star formation. The integrated 
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Fig. 4. Density profiles of satellites, shown in bins of 4" f'- 35 kpc at z ~ 1.8) for clarity. Clockwise from the top left: B- and A",-band flux, 
SFR and stellar mass. The best-fit /3-model, shown with its associated uncertainties by the red line and orange envelope, was derived from the 
unrebinned data. The dashed vertical line marks the average putative virial radius of the halos and the level of contamination by interlopers (i.e., 
“background”) is shown in green. 


stellar mass of UVJ-selected passive and star-forming satellites 
were then individually corrected using iTomczak et al.l (12014h 
mass functions (MF) extrapolated to 10® Mq. A correction 
to the SFR was similarly estimated, based on the main se¬ 
quence parametrization shown above. The corrected total stellar 
mass and SFR are then M^*jot,sat - (9.2 + 2) x 10'° Mq and 
SFR,o,,ja, = (79 + 15) Mq yr“', or respectively ~ 68% and ~ 35% 
of the contribution of the central galaxy. Here we have used 
a “ca nonical” value of 0.8 f or the slope of the main sequence 
(e.g., [Rodighiero et al.ir2014t) . Some rece nt works, on the other 
hand, tend to favor a v alue of near-unity (lAbramson et al.ll2014t 
ISchreiber et al.l I^L5h . If we assume this value, the corrected 
SFR value becomes ~ 54 + 10 Mq y r~*, or ~ 24% of the 
central’s. We note that using lllbert et akl (12013h MFs, derived on 
the same field but from a shallower sample, yields very similar 
values. 


3.1. Far-infrared stacks 

On the other hand, we have extrapolated the MFs to a mass 
range where they are not constrained, and for an environment 
with higher density than the sample from which they were 
defined. Similarly, the slope of the main sequence at this redshift 
is mostly unknown below ~ 10® Mq. Furthermore, extinction- 
corrected SED models might still underestimate SFRs in the 
case of high dust obscuration. The estimated total mass and SFR 
of satellites are therefore uncertain. As an independent check. 


we thus also estimated the total SFR of satellites from available 
FIR maps. Since the resolution of these data (FWHM~ 6 - 20") 
is similar to (or even larger than) the characteristic size of the 
satellite profile, as seen in Fig. |4] an annulus-based analysis 
would likely assign a significant fraction of the IR flux emitted 
by satellites to the central. In addition, satellites close to the 
mass limit are unlikely to be detected in the relatively shallow 
Herschel maps, regardless of their separation from the central. 
To derive the total contribution of satellites to the infrared 
flux of the halos, we have instead stacked, for each band, 
80"x80" cutout images around each central. These stacked 2D 
images were then decomposed into a central point source (for 
the central), a PSF-convolved yS-model centered at the same 
position (for the satellites) and a constant background term. 
For this fit, we have used MIPS 24 fim PSF imag es based on 
obser vations of the GOODS-North held (as used in lElbaz et al.l 
1201 Ih and Herschel PSEs provided by the PEP and HerMES 
collaborations. The parameters of the j8-model were fixed to 
those derived from the catalog-based SER stack and flux uncer¬ 
tainties were estimated through Monte Carlo simulations based 
on parameter errors yielded by the fit, after renormalization so 
that the be at least one per degree of freedom. Eig. |5] shows 
the stacked images in the MIPS 24 jim, HerschelfPKCS 100 and 
150 fim, and Herschel/SPIRE 250 and 350 /im bands, along 
with the best-fit 2D models and residuals. The decomposition 
fails in the case of the SPIRE data, yielding only upper limits 
for the first two bands and providing no meaningful constraint 
to the 500 fim flux of satellites. This is likely due to the high 


Article number, page 5 ofll2l 

































A&A proofs: manuscript no. ms 



Fig. 5. Stacked cutouts (top), best-fit 2D models (middle), and residual 
images (bottom) from the 24/rm Spitzer/MlPS, 100 and 160yum Her- 
schelfPACS, 250 and 350yum Herschel/SPIRE data (left to right; the 
500yum SPIRE image is not shown). The white bar in each top cutout 
has a length of 30". We note that in both the cutouts and model images 
the flux is dominated by the central point-source. 


confusion limit of the instrument, precluding a straightforward 
determination of the background (as shown by the residual 
images in Fig. |5]l, and to the size of the beam being comparable 
to or larger than that of the halos themselves. 


The total infrared luminosity Lik was then derived from the 
resulting SEDs using iMagdis et al.l (l2012h templates convolved 
with the redshift distribution of centrals. We have considered 
both main-sequence and starburst templates, and found that the 
latter perform significantly wors e, as shown in Fig . |6] Convert¬ 
ing L/« into SFR assuming the iKennicutl (Il998h relation, we 
find SFR,o,JR - 176 +11 and 47 + 6 Mq yr for the centrals 
and satellites respectively, consistent with the values derived 
from extinction- and completeness-corrected UV SFRs (which 
we then use in Section H. We note that the total contribution 
of satellites is not sufficient to alter the apparent star formation 
mode (i.e., main-sequence or starburst) of the centrals as 
determined from the HerschelfPACS data. Finally, we can add 
to this value the SFR derived from the uncorrected rest-frame 
UV. Using the B-band flux as a measure of the rest-frame 1500A 
emission, this yields 14 + 3 and 8+1 Mq yr“' for the satellites 
and centrals, respectively. The total SFR of satellites estimated 
from FIR and uncoiTected UV is then 60 + 7 Mq yr"', or ~ 33% 
of the central’s. 


3.2. X-ray observations 

Since our sample is slightly d ifferent and smaller than the one 
used in iBethermin et al.r(l2Q14l) . the average mass of the host ha¬ 
los studied here might be different from that reported in that pa¬ 
per. As previously, we have used deep X-ray observations of the 
COSMO S field by the Chandra and XMM-Newton o bservatories 
(see, e.g.. lFinoguenov et alJl200'^lElvis et alJl2009h to constrain 
the total mass of the halos through a stacking analysis. Of the 
215 centrals in the sample, 14 are directly detected as extended 
sources and 152 are located in zones free from emission. Most of 
the direct detections appear to be consistent with chance associ¬ 
ations along the line of sight with lowe r-redshift galaxy gr oups, 
including 3 that were already known (iGeorge et al.l 1201 lb . We 
have accordingly excluded the “direct detections” from the X- 
ray stack. However, keeping or removing these objects from the 
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Fig. 6. Average far-infrared SED (SpitzerfMlPS 24 pm, HerschelfPACS 
100 and 160 pm, and Herschel/SPIRE 250 and 350 pm) of centrals and 
satellites, with best-fit iMagdis et all (120121) main-sequence templates 
(solid lines). As comparison, best-fit starburst templates are shown as 
dotted lines. Both sets of templates have been broadened according to 
the centrals’ redshift distribution. 


sample has no effect on the rest of the analysis presented in this 
paper and its conclusions. For the sources in regions free from 
detectable emission, we have used the background-subtracted 
and exposure-corrected X-ray image, after subtracting detected 
point sources. The average flux in the 0.5-2 keV band is then 
1.1 X10"'® erg cm"^ s"', detected at 5.3cr. Fig.|2]shows a stacked 
image of these individually undetected objects. For halos in the 
range z ~ 1.5 - 2.5, using the calibrations of iLeauthaud et all 
(1201 oh . this flux corresponds to a rest frame 0.1-2.4 keV lumi¬ 
nosity of 0.8-2.9x10^^ erg s"', an intergalactic medium temper¬ 
ature of ~ 1 keV and a total mass of Mp.nn = 2.1 -2 .4 x 10'^ Mq, 
values similar to those reported in iBethermin et al.l (l2014h . Such 
sources might then soon be indi vidually detectable in de eper X- 
ray surveys such as th e CDF-S (iFinoguenov et al.ll2014h . where 
the applicability of the ILeauthaud et al.l ( 201 Oh scaling relations 
has already been verified for sources with fluxes close to that 
reported here. 


3.3. Comparison with modei predictions 


We have compared the integrated mass and SFR of satellites to 
the predictio ns of different halo occupation distribu t ion (H OD) 
models from Leauthaud et ^(l2012h . lBehroozi et al.l(l2013h . and 
iBethermin et al.l (l2013i hereafter, respectively, LI2, Bhl3 and 
Btl3). The L12 model is shown for its highest defined redshift 
bin (z - 0.74 - 1) while the last two are both evaluated at 
z = 1.8. Fig. |8] shows this comparison for three quantities: the 
ratio of the total stellar mass of satellites to that of the central, 
MsatlMcen, the fraction of stellar mass (central and satellites) to 
total mass, and the ratio of the total SFR of satellites to 

the SFR of the central, SFRja,/SFRce„. We have here used the to¬ 
tal halo masses derived from X-ray stacking, MF-corrected stel¬ 
lar masses and UV+FIR SFRs. All quantities assume the same 
IMF. 


The contribution of satellites to the total stellar mass is 
consistent with the predictions of Btl3 and below L12 at 
z < 1 (this is to be expected, since MsatfMcen decreases with 
increasing redshift in both models). We note that this would 
still be the case if only the measured, uncorrected satellite mass 
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Fig. 7. Stacked X-ray image at the position of the centrals (excluding 
individual detections), from the combined Chandra and XMM-Newton 
0.5 - 2 keV image of the COSMOS field. For reference, the green circle 
has a radius of 15". 



logMh C^s) 


were used. The mass fraction of satellites also appears to be 
compatible with that predicted at z = 0 by rece nt numerical 
simubtions, in the case o f ~ 10'^ Mq halos (iGenel et al.l 

12014 iKravtsov et af] I20l4 . In our probed mass range, the 
stellar mass of the central does not evol ve much with redshift 
dMoster et al.ll2013t iBehroozi et al.ll2()T^ . but it is not obvious 
that this should also be the case for the stellar content of 
satellites. This might therefore suggest that the processes that 
determine the baryon conversion efficiency of the host halo 
also determine to the first order that of the sub-halos. These 
processes can either act within the host halo (what is usually 
thought of when considering environmental effects), or in the 
large-scale structure containing both the central and satellites. 
In the second case, they would then be related to the conformity 
of galaxy prope rties on la rge sca les seen at low redshif t (e.g. , 
iPark et al.ll200^ lAnn et ahlboOSl: iKauffmann et alJl2()in l2013h 
and help synchronize the stellar mass build-up of the central and 
its future satellites, before the satellites merge with the host halo. 

On the other hand, the measured total stellar mass fraction 
is somewhat lower than the predictions ofBtl3 and more 
consistent with Bhl3. This is not entirely surprising, as the for¬ 
mer is optimized to reproduce FIR counts while the latter adopts 
a more sophisticated treatment of the stellar-to-halo mass rela¬ 
tion. The two models also use different stellar mass functions. 
Furthermore, Fig. [8] does not include the systematic uncertainty 
on stellar mass estimates (~ 0.2 - 0.3 dex; see also Section|2]). If 
we take it into account, both the Btl3 and L12 models become 
compatible with the measured value. Finally, the derived SFR 
ratio of satellites and centrals is substantially lower than model 
predictions. This might seem surprising, since the mass ratio is 
itself fully consistent with expectations. On the other hand, the 
total SFR of satellites is, in this model, somewhat sensitive to 
both the behavior of the MF at low masses and the slope of the 
main sequence. For example, if we assume a slope of unity, in¬ 
stead of the value of 0.8 used by Btl3, the predicted SFR ra¬ 
tio would decrease by a factor of ~ 2, making it more consis¬ 
tent with observations. The Btl3 model also adopts a relatively 
simplified treatment of star formation in sub-halos: in the “no¬ 


Fig. 8. Comparison of the integrated properties of hal os in our sample, 
as a fu nction of halo mass, with mode l prediction s from|Bethermi ir et all 
(1201.31 blackl. lLeauthaud et al.lll20T3 . greeni. and lBehroozTeraiTil^oT^ 
red): stellar mass fraction of satellites with respect to the central (top), 
stellar to total mass ratio (middle), and SFR fraction of satellites com¬ 
pared to that of the central (bottom; the dashed line corresponds to a 
case without environmental effects). The black circle and gray diamond 
show, respectively, the MF-corrected and uncorrected values. The latter 
have been shifted to the left for clarity. 

environment” case, the SFR and the quenched fraction are both 
a function of sub-halo mass, while in the other case the model as¬ 
sumes that all satellites of active centrals are themselves active. 
Notably, suppressed (but non-zero) star formation and gradual 
quenching are not considered. 

4. Satellite properties as a function of radius 

In this section, we investigate the variation of the stellar popu¬ 
lation properties of star-forming satellites with distance to the 
central. As in Section |3l we have selected star-forming satellites 
based on t heir rest-frame UVJ c olors, using the high-redshift 
criterion of lWilliams et al.l (l2009l) . Conservatively, we have also 
excluded nominally star-forming objects that are within 0.1 mag 
of the dividing line, so as to avoid possible contamination 
from quiescent satellites. Fig. |3 shows the radial dependency 
of dust extinction, stellar mass, and specific star formation rate 
(SFR/M*, or sSFR). We here look at the variation of median 
values to minimize the effects of outliers. However, because we 
can expect ~ 20% of spurious associations even in the central b 
in (see Fig. O, this measure could still be skewed by interloper 
contamination. To mitigate this, we performed, for each mea¬ 
sured quantity, the following statistical background subtraction: 
in each radial bin within rv,>, we randomly removed a number 
of satellite candidates corresponding to the expected number 
of interlopers, using the background distribution as prior. The 
uncertainties were estimated from the dispersion of median 
values of these background-subtracted distributions. To these 
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values, we have added, as in Section[2 the uncertainties derived 
from bootstrap resampling. This subtraction was performed 
up to the putative virial radius, although the satellite counts 
start becoming consistent with background levels already at 
r > 20"(or ~ 170 kpc; see Fig. [3]). For comparison, Fig.|9]also 
shows the median value prior to the statistical subtraction. 

The stellar mass and extinction of satellites does not vary 
very much with radius, except in the central bins, where the 
median M* and E(B-V) are higher than the background value 
by ~ 0.2 dex and ~ 0.04 mag, respectively. More surprisingly, 
while the SFR density of satellites increases monotonously 
with decreasing halo-centric distance (Fig. S), their median 
sSFR varies significantly within the halo, first exhibiting a 
mild (> 20%, ~ 3cr with respect to the background value) 
rise at 0.5ri,,> (~ 150 kpc), then a more significant decrease 
(~ 40%, 5cr) close to the central (< 50 kpc). Several inter¬ 
pretations, which we discuss below, could account for this effect. 

Normal sSFR variation, as a consequence of the higher 
median mass of satellites at the center since, in the case 
of a non-unity slope for the SFR-M* relation, the sSFR 
is mass-dependent. Fl owever, assuming a slope of ~ 0.8 

(iRodighiero et al.l2014l) and considering the ratio of the sSFR of 
individual galaxies to that of the main sequence at their stellar 
mass does not decrease the significance of the sSFR drop at 
small radii and only slightly that of the excess at 100 - 200 kpc 
(from 3cr to 2cr), as shown in Fig. 0 (bottom right). In fact, a 
very shallow slope of ~ 0.3 would be needed to fully account 
for the observed sSFR decrease. This value seems unlikely for 
~ 10'° Mq galaxie s at z ~ 2 (even in the case of a broken 
power-law relation: IWhitaker et al.l 1201411 and, in the case of a 
single power-law, would reinforce the excess at > 150 kpc. On 
the other hand, a slope value of near-unity would not alter the 
shape of the sSFR variation. 

Stellar population modeling bias or “missed” SFR from 
heavily obscured star formation, e.g., due to a systematic 
underestimation of the extinction-corrected SFR in redder 
galaxies. We have performed a set of simulations using our 
stellar population models with varying extinctions and S/N 
(F(B-V) = 0-1 and S/N> 3, respectively) to test the first 
possibility and quantify the bias to stellar population properties 
in our SED fitting. We find that, when increasing the extinction, 
faint objects will tend to have their stellar mass underestimated 
by ~ 0.05 dex and their reddening and SFR overestimated by 
~ 0.02 mag and < 10%, respectively. These values are within 
the uncertainties of their respective parameters. This is not 
very surprising, as extinction-corrected SFRs derived from 
UV-NIR SEDs have alr eady been found to be quite robust (e.g., 
[Rodighiero et al.ir2014t) and the COSMOS field benefits from a 
large multiwavelength coverage. 

In extreme cases heavily obscured star-forming regions in 
galaxies could be missed entirely by UV-based estimates. In 
this scenario. Fig. |2 could then be interpreted as implying a 
change of star formation mode in satellites as they fall closer 
to the central, toward heavy obscuration. This would result 
in a systematic underestimate of the SFR at small radii. Such 
obscured star formation would however still contribute to the 
integrated rest-frame near- and far-IR light of the galaxies, 
and its influence be detectable in broad-band photometry. As 
shown in Fig. [TOl (left), satellites closer to the central galaxy do 
indeed appear to have slightly redder (~ 0.2 mag) rest-frame 


V-J colors than field galaxies, although the two populations 
still have compatible U-V values and remain within the locus 
of low-extinction star-forming galaxies. This color difference 
could be due to a combination of factors, such as longer 
star formation timescales or higher metallicities (both would 
increase (V-J)o but decrease {U-V)o) in association with higher 
ages or extinction (which increase both (V-J)o and {U-V)o). 
On the other hand, the excellent agreement of the FIR SED of 
satellites with main-sequence models, and between the EIR and 
SED-derived total SERs, suggests that “hidden” s tar formation 
is no t present in significant quantities (see also IZanella et al.l 
1201 5l for similar conclusions). 

Environmental effect on the activity of satellites, from 
interaction with the halo and/or the central. In this case, the 
observed sSER decrease could originate from two different 
galaxy populations: systems with non-zero but suppressed star 
formation, and galaxies where it has recently ceased altogether. 
As the rest-frame UV bands we have used to select and charac¬ 
terize star formation directly trace the light of massive stars, our 
estimates are only sensitive to timescales in excess of 100 Myr 
(e.g., a galaxy can be expected to stay in the star-forming 
locus of the UVJ plane for ~ 300 - 500 Myr after cessation of 
star formation). However, recently quenched systems can be 
easily distinguished from their still active counterparts due to 
the aging of their stellar population affecting bluer bands first. 
Eig. [To] (right), for example, shows two rest-frame UV colors 
of satellites compared to young stellar population models with 
and without ongoing star formation. At z ~ 1.8 the B, V, and r 
filters sample the 1500 - 2300A rest-frame and are thus very 
sensitive to UV light from short-lived massive stars. We see 
no correlation between the shape of the UV continuum and 
distance from the central, with satellites at all radii being on 
average consistent with ongoing star formation. 

We therefore conclude that the observed sSER decrease in 
UVJ-selected star-forming satellites reflects an actual depres¬ 
sion of star formation induced by the group environment. We 
can estimate a lower limit on the timescale of this effect in the 
following way: in a pure free-fall case, a galaxy along a radial 
orbit starting at rv,,> would reach the center of the halo after ~ 
600 - 700 Myr (depending on the concentration of the total mat¬ 
ter distribution, if we assume a NEW profile), reaching velocities 
of ~ 1000 - 2000 km s * and needing only ~ 100 - 150 Myr to 
cross the last 150 kpc, i.e., the radius corresponding to the ob¬ 
served decrease of sSER. If we assume, to the first order, that the 
sSER drop is due to the absence of gas accretion from the satel¬ 
lites’ reservoir s, and that recycling plays a negligible role, we 
find (following lErbll206^ that the time required for the sSER to 
decrease to the observed level would indeed be < 150 Myr. This 
is illustrated in Fig. [TT| where we plot, as a function of radius, 
the diminution of sSFR, assuming the satellites experience no 
gas infall at r < 150 kpc. A more circular orbit would increase 
the time spent by the satellites interacting with the inner halo, 
while a more gradually diminishing gas supply (as well as some 
recycling) would also increase the quenching timescale. 

Several mechanisms can affect the gas reservoirs of galaxies 
in dense environments and induce a diminution of star formation 
(for a review, see, e.g., iBoselli & Gavazzill2006t IPark & Hwand 
|2009|) . Interactions between satellites are here likely not a 
significant driver of galaxy evolution, as the galaxy density 
around individual centrals is relatively low. The minimum 
separation of satellites is ~ 70 - 90 kpc in projection closest 
to the central (where the signal is dominated by real satellites 
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Fig. 9. Median background-subtracted stellar population parameters of UVJ-selected star-forming satellites (open symbols), as a function of 
distance to the central, in bins of 4": stellar mass (top left), extinction (top right), sSFR (bottom left), and sSFR as a fraction of the main- 
sequence value (assuming a slope of 0.8). Median values for the total distribution of star-forming satellites (i.e., without subtracting the background 
distribution) are shown in gray, slightly shifted to the left for convenience. Background levels and related uncertainties, estimated at r > 50", are 
shown in green and the virial radius is indicated by a dashed line. 



Fig. 10. Left. Median rest-frame U-V and V-J colors, derived from SED modeling, of star-forming satellites (filled circles, colored as a function 
of distance to the central) and centrals (open circle and gray density histogram). The purple lines show the expected colors of stellar populations 
models of varying age following a “main-sequence” SFH, assuming E(B-V) < 0.3. Right, median observed B-V and V-r colors of satellites, as a 
function of distance to the central, compared to the expected colors of a star-forming model (main-sequence SFH, black) and a constant SFR one 
observed > 100 Myr after the quenching of star formation (gray). We here include dust extinction with E(B-V) -0-2. The ellipses show the 
dispersion of B-V and V-r values in each radial bin. 
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rather than interlopers; see Fig.O, an order of magnitude larger 
than the typical galaxy s ize in this redshift and mass range (e.g., 
Ivan der Wei et al.ll2014t) . and alr eady above the sc ale at which 
galaxy “harassment” is effective dMoore et al.lll996ll . Because of 
the redshift uncertainties for individual satellites, we can expect 
that the actual distance between them be significantly higher. 
On the other hand, interaction with the hot diffuse intra-halo 
gas, whose presence is confirmed by X-ray stacking, constitutes 
a more plausible source of environmental forcing. The hot 
gas medium can efficiently shut down star formation, mostly 
through hydrodynamical interaction, by either preventing 
further accretion o f cold gas onto the g alaxies (e.g., “starvation”, 
iLarson et al.lll98Cil: iBekki et aDl200:^ or through outright strip¬ 
ping o f the galaxies’ interstellar gas JOunn & Got3ll972l:lNulsenl 
Il982h . These mechanisms are commonly invoked to explain 
general properties of galaxy populations in clusters, such as 
systematic sSFR difference s with respect t o field galaxies (e.g., 
Ivon der Linden et aP 120101: lAlberts k al.l 12014l) and the lack 
thereof. In particular, in massive, high-redshift clusters the sSFR 
of star-forming galaxies does not a ppear to be much cor related 
with cluster-centric distance (e.g., iMuzzin et al.l l2012h . This, 
together wit h the phase-space d istribution of different galaxy 
populations (iMuzzin et al.ll201^ , is viewed as a sign that the 
quenching of star formation in dense environments happens 
on short timescales. The systems studied here probe not only 
a somewhat higher redshift range than the aforementioned 
studies, comparable in fact to the current limit of massive cluster 
samples, but also a mass range that is an order of magnitude 
lower. They are dynamically simpler than large clusters and 
with lower velocities, gas temperatures, and densities. The 
interactions of satellites with their environment should then be 
less violent. Longer interaction timescales might thus explain 
the apparent discrepancy between our analysis, which finds a 
clear sSFR trend, and cluster studies, where such an effect is not 
seen. On the other hand, in the limit case described above (radial 
orbit, no gas infall at < 150 kpc), a galaxy falling toward the 
halo center would have its sSFR decrease by 1 dex in ~ 1 Gyr, 
corresponding to an e-folding time of ~ 0.3 Gyr. This short 
timescale is similar to that inferred for massive clusters and 
consi stent with a fast quenching (see also, e.g., IWetzel et al.l 
l2013h . While constraining the actual mechanisms acting on the 
satellites is beyond the scope of this paper and of the data, we 
note that such timescale is still consistent with either classical 
“starvation” (i.e., mechanical stripping of the gas reservoir; 
iBekki et al.1 [2002h or shock heating of the gas, as predicted 
by hydrodynamical simulations (here, the interaction between 
satellites and their host halo would happen at z < 2.5 in all cases, 
at an epoch when halos of > 10^^ M© are expected to be hot 
and th us prevent efficient c ooling of the g asilDekel & BrrnboirnI 
I 2 OO 6 I) . On the other hand, IZiparo et al.l (1201 3l) report no such 
sSFR gradient in lower redshift groups of similar mass. This 
might reflect a difference between 10^^ M© halos at z ~ 2 
and z ~ 1, in timescales for environmental processes or of 
baryon content. We note however that their highest redshift bin 
(1.2 < z < 1.7) shows a hint of a ~ 0.4 dex sSFR drop similar 
to what we report here, although it is not significant enough 
due to the bin co ntaining only one object of uncertain nature 
dKurk et alJl2009ll . 

On the other hand, the processes described above cannot 
account for the observed sSFR excess in satellites at ~ 150 kpc. 
We can discount galaxy-galaxy interactions for the same rea¬ 
sons, the median minimum separation of satellites in individual 
halos being even higher at > 100 kpc and indistinguishable 


c 
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Fig. 11. Decrease in sSFR in a simple no-infall, no-recycling case as a 
function of distance to the center of the halo (gray, right axis; arbitrar¬ 
ily starting at r = 150 kpc) and tidal perturbation parameter Pg^ (red 
to blue, left axis), for a low-mass satellite. The horizontal line shows 
the threshold value at which the gravita tional influence of the h alo can 
trigger inflow within satellite galaxies dBvrd & Valtonenlll990h . Both 
quantities are shown for NFW concentration parameters ranging from 
c = 1 to c = 40. 


from field levels. Tidal interaction with the halo, however, by 
perturbing the gas already in the galaxies, could accelerate 
star for mation in satellites a nd contribute to clearing their 
gas (see IValentino et al.l 120151 at similar redshift). Following 
iBvrd & ValtonenI (Il99(]l) . we estimate the tidal perturbation 
parameter Pgc - (M/,(r)/M)/(r/A)^, where M, A, r and Mh(r) 
are, respectively, the mass and radius of the satellite, its distance 
from the halo center and the halo mass enclosed within r. For 
Mfiir) we assume a NFW profile with a varying concentration in 
the range c - 1^0 and a total mass given by the X-ray estimate. 
We find that in our case the tidal perturbation starts becoming 
significant {Pgc > 0.006, assuming no stabilizing stellar halo) 
at a distance of r < 200 kpc from the halo center, assuming a 
characteristic sate llite mass of M = 7 x 10 ^ M© (see Fig. |2]i and 
size of A ~ 3 kpc (Ivan der Wei et al.|l2014h . 


4. 1. Passive fraction 

Finally, we note that the observed sSFR decrease at small radii is 
not mirrored by an increase in the number of quiescent satellites 
near the central. In Fig. [12] we show the ratio of UVJ-selected 
quiescent satellites to the total number of satellites in each radial 
bin. We performed the same statistical background subtraction 
as described above, using the color distributions of the satellites 
as priors and estimating the quiescent fraction for each random 
trial. The fraction of quiescent galaxies is close to 20% at large 
radii and appears to decrease slightly at r < 20". This valu e and 
trend are similar to those derived by iHartlev et al.l (1201 5h from 
a slightly lower redshift sample. On the other hand, if we adopt 
a slightly more stringent criterion, by adding a 0.1 mag mar¬ 
gin (see Fig. ITOl l and selecting only the redder UVJ-quiescent 
galaxies, the background quiescent fraction drops to ~10% and 
the trend at small radii disappears. This suggests that, at least 
in our case, it is mostly due to objects close to, or straddling, 
the dividing line between the passive and star-forming loci. The 
stellar mass distributions of both quiescent samples are not sig¬ 
nificantly different, however. The absence of a clear number ex- 
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Fig. 12. Fraction of quiescent satellites as a function of radius. The 
satellites were selected according to the UVJ criterion (open symbols), 
and adding a margin of 0.1 mag in both (U-V)o and (V'-y)o (gray sym¬ 
bols). Values within the virial radius, as shown by the dashed vertical 
line, were derived from background-subtracted distributions using the 
same procedure as in Fig. [O] Background levels and related uncertain¬ 
ties, estimated at r > 50", are shown in green. 
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cess could seem counterintuitive, considering the sSFR variation 
described in Section 01 On the other hand, the appearance of an 
obvious quiescent galaxy population takes time. For it to hap¬ 
pen in this case, the quenching of star formation in the central 
satellites would have to have started at z ~ 2.5, if we assume 
a time span of 1 Gyr for a 1 dex sSFR decrease as discussed 
above. This would in turn imply that the environmental condi¬ 
tions responsible for it (e.g., a hot halo) be already in place at 
this epoch. We can infer that this was not the case in the type of 
halos investigated here. 


5. Conclusions 

Low-mass structures traced by massive galaxies, while more dif¬ 
ficult to confirm individually, can be efficiently selected statisti¬ 
cally. At high redshift, they can offer a more accessible window 
to galaxy evolution in dense environments than galaxy clusters, 
their high abundance compensating for the lower galaxy num¬ 
ber density and environmental bias, even in a relatively limited 
area. We have here taken advantage of the wealth and depth of 
photometric data available on the COSMOS field to study the 
distribution and properties of star-forming satellites associated 
with massive galaxies on the main sequence of star formation, 
as tracers of group-size halos of mass ~ 2 - 3 x 10^^ M©. We 
have constructed a sample of massive star-forming galaxies at 
(z) = 1.8, selecting only objects without close neighbors of com¬ 
parable mass so that they be putatively central to their host halo. 
We have verified the average total mass of said halos thanks to 
deep Chandra and XMM data, and found it to be < 3 x 10^^ M©. 
Using the recently released matched photometric catalogs for the 
COSMOS field, we have derived stellar population parameters 
for both centrals and satellites. Our conclusions are the follow¬ 
ing; 


- we have estimated the contribution of satellite galaxies to the 
stellar mass and SFR of the systems at, respectively, ~ 68% 
and ~ 25 - 35% of the stellar mass and SFR of the central 
galaxy (or ~ 40% and < 25% of the total stellar mass and 


SFR), after correcting for the completeness limit of the sam¬ 
ple. The stellar mass fraction of satellites with respect to the 
central is found to be consistent with the predictions of HOD 
models, as is the total stellar mass to halo mass ratio. On the 
other hand, the observed total SFR of satellites appears to be 
a factor of ~ 2 - 3 lower than model predictions. This might 
be related to the relatively simple treatment of star formation 
in sub-halos adopted by our chosen model, or to assumptions 
on the behavior of the main sequence of star formation at low 
stellar mass. 

- we have also independently estimated the SFR of satellites 
and centrals from stacked FIR data, by separating their con¬ 
tributions through source decomposition. The SED thus de¬ 
rived is well-fitted by a main-sequence template and yields a 
SFR of ~ 47 M© yr“', consistent with the UV-NIR estimate. 
This also suggests an absence of significant heavily obscured 
star formation (e.g., starbursts) in the satellite population. 

- finally, we have probed the radial dependence of the prop¬ 
erties of star-forming satellites. We find significant variation 
of their sSFR within the virial radius, with a marginal excess 
at r ~ 150 kpc followed by sharper drop at r < 100 kpc. 
This suggests that the group environment acts differently on 
star-forming galaxies within rv,,> depending on their distance 
to the center, enhancing star formation slightly at larger radii 
while quenching it with a timescale of > 300 Myr closer 
to the center. In the first order, this is consistent with desta¬ 
bilization of galactic gas by the halo potential followed by 
prevention of further gas accretion, as the galaxy falls closer 
to the center of the halo. 

On the other hand, the use of photometric data not only 
implies some amount of back- and foreground contamination, 
but also precludes knowledge of important quantities, such as 
the instantaneous star formation rate and metallicity, that would 
more precisely constrain the mechanisms of galaxy evolution 
in these halos. Wide-held, high-coverage spectroscopic instru¬ 
ments (e.g., large integral held units such as MUSE) and, later, 
“all-in-one” large-scale surveys (e.g., Euclid and WEIRST), 
should allow for a dramatic improvement in statistics and red- 
shift resolution, especially around the critical epoch of galaxy 
and cluster progenitor build-up at z ~ 2. 
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